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ABSTRACT 

We present the first resolved study of the radio continuum properties of 
IZwl8, the dwarf galaxy with the lowest known nebular metal abundance in 
the local universe. New Very Large Array radio continuum images at 20 and 3.6 
cm are compared to various Hubble Space Telescope images, and we find a strik- 
ing morphological similarity between high resolution Ha and short wavelength 
radio continuum emission, especially in the Ha shell in the northwest region. We 
separate thermal and nonthermal components of the emission, and find a large 
synchrotron halo surrounding the galaxy. Comparison between Ha and X-band 
fluxes suggests that the emission at 3.6 cm is dominated by thermal processes; 
an additional synchrotron component dominates the flux at 20 cm and produces 
a modest fraction of the detected flux at 3.6 cm. The fluxes of three of the four 
major emission peaks show a mix of thermal and nonthermal processes, while 
one shows a nearly flat spectral index. The strong synchrotron component ar- 
gues for active star formation throughout the disk for at least the last ~ 30 Myr. 
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These sensitive observations provide a new, detailed view of the nature of radio 
continuum emission in the very low metallicity interstellar medium. Comparing 
with the literature, the role of metallicity in the evolution of radio continuum 
emission seems to be secondary to other factors such as recent star formation 
history and the presence or absence of outflows from star formation regions. 

Subject headings: galaxies: evolution — galaxies: dwarf — galaxies: individual 
(I Zw 18) — radio continuum: galaxies 

1. Introduction 

With the lowest nebular metallicity known in the local universe [(O/H) ~ 0.02 (O/H) ; 
Skillman & Kennicutt 1993], IZwl8 plays a key role in the study of galaxy evolution. Un- 
derstanding the nature of the star-forming fragments predicted at high redshift requires a 
detailed study of nearby potentially "young" galaxies. Estimates of the age of the stellar 
population in I Zw 18 based on single-star photometry suggest an age of less than 1 Gyr (e.g., 
Hunter & Thronson 1995; Aloisi etal. 1999; Ostlin 2000; Izotov & Thuan 2004). 

An alternative probe of the unobscured, recent (< 30 Myr) star formation in a galaxy 
can be obtained with radio continuum observations. Three types of radio continuum emis- 
sion are seen in local star-forming galaxies (see Condon 1992 for a detailed review). First, 
thermal free-free emission (characterized by a nearly flat spectral index, a ~ —0.1, where 
~ u a ) is proportional to the total Lyman continuum flux and hence can be used as an 
accurate star formation rate estimator. Second, nonthermal synchrotron emission, produced 
in supernovae (SNe) explosions and remnants, shows a more negative spectral index (—1.2 
< a < —0.4). Finally, free-free absorption can occur in young, dense, heavily-embedded 
clusters that show an inverted spectrum (a > 0.0; Turner etal. 1998, Kobulnicky & Johnson 
1999, Johnson & Kobulnicky 2003). Distinguishing between these emission mechanisms al- 
lows one to constrain the ages of the emitting components: inverted sources are very young 
(ages ~ 10 5 - 10 6 yr; Johnson etal. 2001); thermal regions have ages comparable to those of 
H II regions (^ 10 Myr); synchrotron radiation dominates in regions older than this. 

Considering the amount of observational attention that I Zw 18 has received, it is some- 
what surprising that no detailed investigation of its radio continuum properties exists in the 
literature. There have been two previously published radio continuum detections, but both 
have been extracted from spectral line observations (Lequeux & Viallefond 1980; van Zee 
etal. 1998). In order to study the nature of radio continuum emission at the lowest avail- 
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able nebular metallicity, we have obtained sensitive Very Large Array (VLA) 1 observations 
of IZwl8. We compare these data to Hubble Space Telescope (HST) 2 images of the stellar 
and nebular emission. Throughout this paper, we adopt a distance of 12.6 Mpc for IZwl8 
(Ostlin 2000); however, most of our results (spectral indices, comparison to HST imaging) 
do not require a tightly-constrained distance estimate. 

2. Observations and Data Reduction 

2.1. VLA Radio Continuum Data 

VLA radio continuum imaging was obtained using the A, B, and C arrays on 2003, 
August 3, and 2004, January 20 and April 20 (total on-source integration time = 15.1 hours) 
for program AC681. All reductions were performed using the Astronomical Image Processing 
System (AIPS) package. First, interference and bad data were removed. Flux, gain and 
phase calibrations were then applied, derived from observations of 1331+305 (3c286; primary 
calibrator) and 0921+622 (secondary calibrator). The L-band (20 cm) and X-band (3.6 cm) 
observations were obtained in two arrays, and these u-v databases were concatenated. The 
calibrated u-v data were then imaged and cleaned to produce matched-beam images at each 
frequency that balanced resolution and sensitivity. We analyze a set of "high" resolution (L- 
band: beam = 2.18" x 2.01", rms = 11.6 /iJybeam -1 ; X-band: beam = 2.17" x 2.07", rms = 
7.4 ^Jybeam -1 ) and a set of "low" resolution (beam = 4.5", L-band rms = 19 iJybeam -1 , 
X-band rms = 12 /zJybeam -1 ) images. 

2.2. HST Imaging 

We compare our radio observations to archival HST observations of I Zw 18 from pro- 
grams GO-5434 (P.I. Dufour), GO-6536 (P.I. Skillman), and GO-9400 (P.I. Thuan). The 
detailed handling of the data from programs 6536 and 5434 are presented in Cannon et al. 
(2002). We use the same narrowband Ha image presented in that work for analysis here 
[total Ha flux = (3.26 + 0.3) x 10~ 13 ergs -1 cm -2 , implying a current star formation rate of 



1 The National Radio Astronomy Observatory is a facility of the National Science Foundation operated 
under cooperative agreement by Associated Universities, Inc. 

2 Bascd on observations with the NASA/ESA Hubble Space Telescope, obtained at the Space Telescope 
Science Institute, which is operated by the Association of Universities for Research in Astronomy Inc. under 
NASA contract No. NAS5-26555. 
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~ 0.05 Moyr^ 1 (Kennicutt etal. 1994)]. From program 9400, we use an HST/ ACS F555W 
image for comparison with the stellar continuum. Based on the positions of HST Guide Stars 
in the ACS field and the agreement of the radio and Ha peaks, we conservatively estimate 
the uncertainty in the coordinate solution to be better than 0.5" rms. 



3. The Nature of Radio Continuum Emission 

3.1. Global Flux Measurements 

In Figure 1 we present the low-resolution images overlaid with contours at various levels. 
We measure global flux densities (by integrating continuous emission surrounding the galaxy 
after blanking at the 3 a level) at L- and X-bands of 1.79 ±0.18 mJy and 0.78 ±0.08 mJy, 
respectively, corresponding to a global spectral index of —0.46 ±0.06 (see Table 1). The 
continuum emission is more spatially extended at 20 cm than at 3.6 cm (although the latter 
data are significantly more sensitive); indeed, emission above the 3a level is detected out to 
~ 12" (deconvolved) northwest of the emission peak at L-band (corresponding to a linear 
distance of ~ 730 pc at the adopted distance), but only to ~ 7.5" at X-band. As argued in 
the next section, this is the result of a synchrotron halo surrounding the system. 



3.2. Thermal vs. Nonthermal Decomposition 

Since the thermal radio continuum luminosity is proportional to the total photoioniza- 
tion rate, there exists a well-defined relation between radio continuum luminosity and other 
direct star formation rate indicators (e.g., Ha emission). In systems with little extinction in 
the optical, this relation can be used to estimate the expected thermal fraction of detected 
radio continuum emission. In IZwl8, the low detected extinctions (Ay < 0.5 mag.; see 
Cannon et al. 2002) imply that Ha provides a nearly complete census of the ongoing star 
formation. 

We can estimate the expected thermal components at each frequency using the relations 
presented in Caplan & Deharveng (1986). For a purely thermal radio continuum source in 
the absence of extinction, the Ha flux and the flux density at a given frequency are related 
via: 

j Ha _ (8.67xl0- 9 )(T)-°- 44 

j v (10.811 + 1.5 • ln(T) - ln(v)) U 



- 5 - 



where ]ua is the flux at Ha in ergs -1 cm -2 , ] u is the flux density at the given frequency in Jy, 
T is the electron temperature in units of 10 4 K, and v is the observed frequency in GHz. We 
adopt an intermediate electron temperature from Skillman & Kennicutt (1993) of 18500 K 
for the entire galaxy. Using the global Ha flux measured from the HST image (see Table 1), 
we expect flux densities of 0.47 mJy at X-band and 0.56 mJy at L-band. Comparison with 
the measured X-band value of 0.78 mJy suggests that more than half of the 3.6 cm emission 
is thermal; conversely, only ~ 30% of the emission at L-band is expected to be thermal. 

Given this expected value for thermal emission at 20 cm, > 70% of the detected emis- 
sion can be attributed to a synchrotron component. If this emission has a characteristic 
synchrotron spectral index (a = —0.8), then this implies a synchrotron contribution to the 
detected flux at X-band of 0.29 mJy. This value can be compared with the difference be- 
tween total and (predicted) thermal fluxes, 0.31 mJy; the agreement suggests that such a 
synchrotron component fits the data quite well. 

3.3. Properties of Individual Radio Emission Peaks 

We present higher resolution images of the radio continuum emission from I Zw 18 in 
Figures 2 and 3, as well as overlays of the emission on HST Ha and V-band images. We 
resolve the emission into 4 distinct peaks at this resolution. Examining these images, there is 
a striking correlation between Ha and radio continuum emission at this sensitivity level; in all 
areas of high surface brightness Ha emission, strong radio continuum emission is prevalent. 
Note also that the radio continuum emission closely follows the diffuse Ha shell morphology 
in the northwest region. 

The total flux densities of the individual peaks were measured using identical apertures 
the size of the beam or larger, centered on the peaks in the X-band image (see Figure 2 and 
Table 1). The individual radio continuum peaks are found to have a mix of thermal and 
nonthermal spectral indices. Three of the sources (I Zw 18 SE, NW-A, NW-C) have spectral 
indices a ~ —0.27, implying a mix of thermal and nonthermal components. IZwl8 NW-B, 
however, is consistent with nearly purely thermal emission (a = —0.04 ±0.06). 

Sources I Zw 18 SE and NW-A are clearly associated with collections of the high surface 
brightness H II regions of the galaxy. The high (aperture integrated) Ha equivalent widths 
found in Cannon et al. (2002) imply strong active star formation, while their nonthermal radio 
continua suggest that the areas are also rich in massive star feedback via SNe explosions and 
remnants. This suggests that the star formation intensity was strong ~ 20-30 Myr ago, as 
well as at the present epoch. This is in good agreement with the results of resolved stellar 
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population analysis, where recent investigations have found a star formation rate that was 
elevated during the last 10-100 Myr compared to that at present (Aloisi etal. 1999; Ostlin 
2000; Izotov & Thuan 2004). 

4. Discussion and Conclusions 

The nature of radio continuum emission at very low metallicities is only beginning to be 
explored. Comparing our results with those on the second-most metal poor galaxy known, 
SBS 0335—052 (Hunt etal. 2004), we find that both systems have a mix of nonthermal and 
thermal components, indicative of active massive star formation over the last ~ 30 Myr. 
However, the radio spectrum of SBS 0335—052 has a thermal component that manifests 
itself in considerable absorption at L-band and a radio turnover between 1.4 and 4.8 GHz. 
While the present dataset has fewer frequencies, we would be sensitive to a similar spectral 
turnover in IZwl8 (as this would result in positive, rather than negative, spectral indices 
between 3.6 and 20cm). 

The differences in the radio continuum properties of I Zw 18 and SBS 0335—052 are most 
easily explained as the results of the more vigorous star formation event that is underway 
in the latter. IZwl8 appears to be undergoing star formation in discrete stellar clusters 
and H II regions [see Table 1 for Ha luminosities of the radio emission peaks studied here, 
and Cannon etal. (2002) for those of other clusters], and infrared imaging (e.g., Ostlin 2000; 
Ostlin & Mouchine 2005) does not suggest embedded regions that would be able to give rise 
to "super star clusters" (SSCs; see O'Connell 2004 for a recent review of SSC properties). 
SBS 0335—052, on the other hand, hosts multiple SSCs (Thuan etal. 1997) and regions 
of exceptionally high extinction and embedded star formation for a very metal poor galaxy 
(Hunt etal. 2001; Plante & Sauvage 2002). Given that complex (and often inverted) spectral 
indices typically arise from such regions (Turner etal. 1998; Kobulnicky & Johnson 1999; 
Johnson & Kobulnicky 2003), it is perhaps not surprising that the radio spectra of IZwl8 
and SBS 0335-052 appear to differ. 

To summarize, we have detected thermal and nonthermal emission in I Zw 18, including 
a synchrotron halo. This low-metallicity, low-dust system shows a radio continuum spec- 
trum between 20 and 3.6 cm that is typical of star-forming dwarf galaxies; however, the 
radio properties differ markedly from those found in the second most metal-poor galaxy, 
SBS 0335—052. Although the sample is small, the variety of radio continua found at very 
low metallicities suggests that there is not a strong metallicity dependence for the produc- 
tion of radio continuum emission. Indeed, it is not necessary that one should exist; thermal 
continua are produced by stars that have formed, while nonthermal continua are typically 
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the result of the evolution of massive stars. While there certainly exist differences in the 
way that the stars form at different metallicities, once the stars have formed, radio contin- 
uum emission appears to be ubiquitous (though complex and dependent on other factors) in 
massive star formation regions regardless of metallicity. 

The authors would like to thank Chip Kobulnicky for helpful discussions during the 
planning stages of this project, and the anonymous referee for helpful comments. 
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Fig. 1. — Low-resolution X-band (a) and L-Band (b) images of IZwl8, overlaid with con- 
tours at the (—4, —2, 2, 4, 8, 16) a levels (X-band 1 a = 12 /iJybeam -1 ; L-band 1 a = 
19 yuJybeam" 1 ). The beam size is 4.5" x 4.5" and is shown at the bottom left. The black 
box in (b) outlines the approximate field of view shown in Figures 2 and 3. Note that the 
diffuse emission is more spatially extended at 20 cm than at 3.6 cm. 




Fig. 2. — High-resolution X-Band emission contours at the —3, 3, 4.5, 6, 7.5, 9, 10.5, 12, 
& 13.5ct levels (la = 7.4 //Jybeam" 1 ), superposed on a grey-scale image of the emission 
(a), on an HST/ ACS F555W (V) image (b), and on an HST/WFPC2 F656N (continuum- 
subtracted Ha) image (c). Labels and the white circles in (a) denote the apertures used to 
measure fluxes of individual peaks in the radio continuum and in Ha; see § 3 and Table 1. 
The beam size is 2.17" x 2.07", and is shown at upper left. The field of view (with north up 
and east to the left) is ~ 0.9 x 0.8 kpc at the adopted distance of 12.6 Mpc. 



Fig. 3. — Same as Figure 2, but for the L-band high-resolution data (la = 11.6 /iJybeam^ 1 ). 
Comparing to the X-band emission at similar resolution, it is clear that the individual emis- 
sion peaks are less well-defined; we interpret this as the result of the additional synchrotron 
component at L-band. The beam size is 2.18" x 2.04", and is shown at upper left. 
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Table 1. Properties of Radio Continuum Emission in IZwl8 



Parameter 


IZwl8 SE a 


I Zw 18 NW-A a 


I Zw 18 NW-B a 


IZwl8 NW-C a 


I Zw 18 Total Galaxy 


R.A. (J2000) a 


09:34:02.36 


09:34:02.10 


09:34:01.85 


09:34:01.82 


09:34:02.0 


Dec. (J2000) a 


+55:14:23.10 


+55:14:28.06 


+55:14:26.00 


+55:14:29.06 


+55:14:28 


Aperture Radius (") 


1.75 


1.5 


1.1 


1.1 




SLBand OJy) 


0.17 + 0.02 


0.20 ±0.02 


0.050 + 0.005 


0.079 + 0.008 


1.79 + 0.18 


SxBand (mjy) 


0.11 + 0.01 


0.12±0.01 


0.047 ±0.005 


0.050 + 0.005 


0.78 + 0.08 


a (S„ ~ v a ) 


-0.27 ±0.06 


-0.29 + 0.06 


-0.04 + 0.06 


-0.26 + 0.06 


-0.46 + 0.06 


Ha Flux (10~ 14 ergs" 1 cm- 2 ) 


7.7 ±0.8 


7.8 ±0.8 


1.3 + 0.2 


2.9 + 0.3 


32.6 + 3.3 



a See Figure 1 for positions. The central position of radio continuum emission peaks are derived from Gaussian fitting. The 
coordinates of the I Zw 18 system are taken from the NASA Extragalactic Database. 



